A locally quadratic Glimm functional and sharp 
convergence rate of the Glimm scheme for 
nonlinear hyperbolic systems 

Fabio Ancona* Andrea Marson^ 
April 19, 2009 

September 1 2008 
Abstract 

Consider the Cauchy problem for a strictly hyperbolic, N x N quasilinear 
system in one space dimension 

Ut + A{u)Uj: = 0, u{0,x) = uix), (1) 

where u k-> A{u) is a smooth matrix-valued map, and the initial data u is 
assumed to have small total variation. We investigate the rate of convergence 
of approximate solutions of (1) constructed by the Glimm scheme, under 
the assumption that, letting Afc(u), rk{u) denote the k-ih eigenvalue and a 
corresponding eigenvector of A{u), respectively, for each fc-th characteristic 
family the linearly degenerate manifold 

Mk^ {uen : VAfe(u) ■ rk{u) = O} 

is either the whole space, or it is empty, or it consists of a finite number 
of smooth, A'^— 1-dimensional, connected, manifolds that are transversal to 
the characteristic vector field rk ■ We introduce a Glimm type functional 
which is the sum of the cubic interaction potential defined in [6], and of a 
quadratic term that takes into account interactions of waves of the same 
family with strength smaller than some fixed threshold parameter. Relying 
on an adapted wave tracing method, and on the decrease amount of such 
a functional, we obtain the same type of error estimates valid for Glimm 
approximate solutions of hyperbolic systems satisfying the classical Lax as- 
sumptions of genuine nonlinearity or linear degeneracy of the characteristic 
families. 
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1 Introduction 

Consider the Cauchy problem for a general system of hyperbolic conservation laws 
in one space dimension 

ut + F{u), = 0, (1.1) 
u{0,x) = u{x). (1.2) 

Here the vector u = u{t,x) — (ui{t, x), . . . ,U]y{t, x)'j represents the conserved 
quantities, while the components of the vector valued function 

Fiu) = {Fi{u),...,FNiu)) 

are the corresponding fluxes. We assume that the flux function is a smooth map 
defined on a domain C R^, and that the system (|l.ip is strictly hyperbolic, i.e. 
that the Jacobian matrix A{u) — DF(u) has N real distinct eigenvalues 

Xi{u) < ■ ■ ■ < Xn{u) Vu. (1.3) 

Denote with ri{u), . . . ,ri^{u) a corresponding basis of right eigenvectors. Hyper- 
bolic equations in conservation form physically arise in several contexts. A primary 
example of such systems is provided by the Euler equations of non- viscous gases, 
see P] . 

It is well known that, because of the nonlinear dependence of the characteristic 
speeds Xk{u) on the state variable u, classical solutions to (|l.ip can develop dis- 
continuities (shock wave) in finite time, no matter of the regularity of the initial 
data. Therefore, in order to construct solutions globally defined in time, one must 
consider weak solutions interpreting the equation (jl.ip in a distributional sense. 
Moreover, for sake of uniqueness, an entropy criterion for admissibility is usually 
added to rule out non-physical discontinuities. In [20l [21] T.P. Liu proposed the 
following admissibility criterion valid for weak solutions to general systems of con- 
servation laws, that generalizes the classical stability condition introduced by Lax 
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Definition 1.1 A shock discontinuity of the k-th family {u^ , u^), traveling with 
speed ak [u^ , u^] , is Liu admissible if, for any state u lying on the Hugoniot curve 
Sk[u^] between and , the shock speed crfc[it^,u] of the discontinuity {u^,u) 
satisfies 

ak[u^,u]>ak[u^,u'^]. (1.4) 

The existence of global weak admissible solutions to (|l.ip - (ll.2p with small total 
variation was first established in the celebrated paper of Glimm |14j under the 
additional assumption that each characteristic field be either linearly degenerate 
(LD), so that 

VXk{u)-rk{u)=^0 Vu, (1.5) 
or else genuinely nonlinear (GNL) i.e. 

VXkiu)-rkiu)^0 Vu. (1.6) 

A random choice method, the Glimm scheme, was introduced in [T7 to construct 
approximate solutions of the general Cauchy problem p.ip - p.2p by piecing to- 
gether solutions of several Riemann problems, i.e. Cauchy problems whose initial 
data are piecewise constant with a single jump at the origin 



u 



^ if X < 



<0,x)^-> „ (1.7) 



u 



if X > 



Using a nonlinear functional introduced by Glimm, that measures the nonlinear 
coupling of waves in the solution, one can establish a-priori bounds on the total 
variation of a family of approximate solutions. These uniform estimates then yield 
the convergence of a sequence of approximate solutions to the weak admissible so- 
lution of p.ip - (|1.2p . The existence theory for the Cauchy problem (|l.ip - (|1.2p based 
on a Glimm scheme was extended by Liu P^, Liu and Yang (23j . and by Iguchy 
and LeFloch [17j to the case of systems with non genuinely nonlinear (NGNL) 
characteristic families whose flux function satisfy the more general assumption: 

(H) The vector valued function F is C'^, and for each k e {1, . . . . iV}-th charac- 
teristic family the linearly degenerate manifold 

Mk^{uen : VAfc(M)Tfe(u)-0} (1.8) 

is either empty (GNL characteristic field), or it is the whole space (LD char- 
acteristic field), or it consists of a finite number < M of smooth, A^— 1- 
dimensional, connected, manifolds, and there holds 

V(VAfcTfc)(w)-rfe(ii)^0 WueMk- (1.9) 



Aim of the present paper is to provide a sharp convergence rate for approximate 
solutions obtained by the Glimm scheme valid for strictly hyperbolic systems of 
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conservation laws satisfying the assumption (H). We recall that in the Glimm 
scheme, one works with a fix grid in the t-x plane, with mesh sizes At, Ax. An 
approximate solution of p.ip - (|1.2p is then constructed as follows. By possibly 
performing a linear change of coordinates in the t-x plane, we may assume that 
the characteristic speeds Afe(it), 1 < fc < A^, take values in the interval [0, 1], for all 
u £ V,. Then, choose At = Ax = e, and let {^^j^gN C [0, 1] be an equidistributed 
sequence of numbers, which thus satisfies the condition 



lim 



n — 1 



£=0 



= 



V A € [0, 1] , 



(1.10) 



where X[o,a] denotes the characteristic function of the interval [0, A]. On the initial 
strip < i < £, is defined as the exact solution of p.ip . with starting condition 

M^(0,x) =u((j+0o)£) Va;e]j£, {] + l)e[. 

Next, assuming that has been constructed for t £ [0, ie[, on the strip ie < t < 
{i + l)e , is defined as the exact solution of (II. ip . with starting condition 

-u^(ie, x) = , (j + 6i)e) V x e ]je, {j + 1)£ [ . 

Relying on uniform a-priori bounds on the total variation, we thus define induc- 
tively the approximate solution u^(i, •) for all t > 0. 

One can repeat this construction with the same values 9i for each time interval 
[ie, (i + l)e[ , and letting the mesh size e tend to zero. Hence, we obtain a sequence 
of approximate solutions which converge, by compactness, to some limit function u 
that is shown to be a weak admissible solution of p.ip - (|1.2p (cfr. [H]). In order to 
derive an accurate estimate of the convergence rate of the approximate solutions, 
it was introduced in [TU] an equidistributed sequence {6'^}fgN C [0, 1] enjoying the 
following property. For any given < m < n, define the discrepancy of the set 
{9m, ■ ■ ■ , dn-i} as 



Dm,n = sup 

Ae[o,i] 



A- 



m<£<n 



Then, there holds 



Drn.n < 0{1) 



1 + log(n — to) 



V n > TO > 1 , 



(1.11) 



(1.12) 



Here, and throughout the paper, 0(1) denotes a uniformly bounded quantity, while 
we will use the Landau symbol o(l) to indicate a quantity that approaches zero as 
e —^ 0. Relying on the existence of a Lipschitz continuous semigroup of solutions 
generated by (jl.ip . compatible with the solutions of the Riemann problems, it 
was proved in fTU] that, for systems with GNL or LD characteristic fields, the 
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convergence rate of the Glimm approximate solutions constructed in connection 
with a sequence enjoying the property (|1.12p is o(l)-y^ | lne|. In the case of general 
systems satisfying the assumption (H), it was derived in [IB] an estimate of the 
norm of the error in the Glimm approximate solutions of the order o(l)-v/e | lne|. 

In the present paper, we improve this result by establishing the same con- 
vergence rate of the approximate solutions generated by the Glimm scheme for 
systems satisfying the assumption (H) as in the case of systems with GNL or LD 
characteristic fields. Namely, our result is the following. 

Theorem 1.2 Let F be a map from a domain fi C into satisfying the 
assumption (H), and assume that the system lil.l]) is strictly hyperbolic. Given an 
initial datum u with small total variation, let u{t, •) be the unique Liu admissible 
solution of il.l\) - n~^} . Let {Ok}k&'i C [0,1] be a sequence satisfying il.lS]) and 
construct the corresponding Glimm approximate solution of U.1] )- [T^) with 
mesh sizes Ax = At — e. Then, for every T > there holds 

^^^mr,.) uiT.)U. (1.13) 
e^o ^/e|log£| 

and the limit is uniform w.r.t. u as long as Tot.Var.(?l) remains uniformly small. 

Our result applies more generally to strictly hyperbolic N x N quasilincar 
systems 

Ut+A{u)u^^O, (1.14) 

not necessarily in conservation form, where ^ is a matrix valued map defined 
from a domain C into M^^^(R), whose eigenvalues A^, k e {1, . . . , A^}, 
satisfy the assumption stated in (H). Indeed, one may alternatively assume that 
A : n ^ M^><^(R) is a C^-i map, and that for each NGNL k € {1,...,N}- 
th characteristic family the linearly degenerate manifold A4k consists of a finite 
number of connected manifolds Aik,h, that are either iV— 1-dimensional as in (H) 
or A''-dimensional with a similar condition to (|1.9p (cfr. Remark 13.61 in § [3] and 
Remark 16.11 in § [5]) . 

In fact, the fundamental paper of Bianchini and Bressan [7] shows that, for 
any C^'^ map A : fl ^ M^^^(R) with strictly hyperbolic values, (jl.lip generates 
a unique (up to the domain) Lipschitz continuous semigroup {St : t > 0} of 
vanishing viscosity solutions obtained as the (unique) limits of solutions to the 
(artificial) viscous parabolic approximation 

Ut+ A{u)Ux^ tJ.Uxx, (l-16)/i 

when the viscosity coefficient /i ^ 0. The trajectories of such a semigroup start- 
ing from piecewise constant initial data locally coincide with the "admissible" 
solution of each Riemann problem determined by the jumps in the initial data. 
Moreover, any limit of Glimm approximations coincides with the corresponding 
trajectory of the semigroup generated by (|1.14p . In particular, in the conserva- 
tive case where A{u) = DF{u) every vanishing viscosity solution of the Cauchy 
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problem (|1.14p - (ll.2p provides a weak solution of p.ip - p.2p satisfying the Liu ad- 
missibility conditions (|1.4p . 

The proof of the error bound (jl.l3|) follows the same strategy adopted in [10] , 
relying on the careful analysis of the structure of the solution for systems satisfying 
the assumption (H), developed by T.P. Liu and T. Yang in [551 Indeed, 
to estimate the distance between a Lipschitz continuous (in time) approximate 
solutions w of (|1.14p and the corresponding exact solution one would like to use 
the error bound [9 

f'^ \\w(t + h) - Shw(t)L ^ 

\\w{T) ^ Stw{0)\l , < L liminf^ , ^ ^"^ dt , (L17) 

" Jo '1^0+ h 

where L denotes a Lipschitz constant of the semigroup S generated by (II. 
However, for approximate solutions constructed by the Glimm scheme, a direct 
application of this formula is of little help because of the additional errors intro- 
duced by the restarting procedures at times ti = ie. For this reason, following the 
wave tracing analysis in [23j , it is useful to partition the elementary waves present 
in the approximate solution, say in a time interval [ti , T2] , into virtual waves that 
can be either traced back from T2 to ri, or are canceled or generated by inter- 
actions occurring in [ti , T2\ ■ Thanks to the simplified wave pattern associated to 
this partition, one can construct a front tracking approximation having the same 
initial and terminal values as the Glimm approximation, and thus establish (|1.13p 
relying on p.l7p . 

As one would expect, the presence of elementary waves with various composite 
wave patterns for systems satisfying the assumption (H), requires a careful analysis 
of the errors introduced by this wave-partition algorithm. As customary, the 
change of wave-size and wave-speeds when an interaction takes place is controlled 
by a Glimm functional that measures the potential interaction of waves in the 
solution. 

For general strictly hyperbolic systems (|1.14p satisfying the assumption (H), 
several nonlinear functionals were introduced in [551 1131 dZl IS]) consisting of a 
standard Glimm quadratic functional, for the interaction of waves of different 
families, and of a cubic functional measuring the potential interaction between 
waves of the same family. This cubic part of the functional is defined in terms 
of the strengths of any pair of waves of the same family and of the absolute 
value of the angle between them [B] (or of the positive part of the angle between 
two waves [551|5^). Such functionals work perfectly to establish uniform a- priori 
bounds on the total variation of the solution, but are not effective to control the 
quadratic order error produced by the change of wave speeds for interactions of 
waves of the same family, of arbitrarily small sizes. 

On the other hand, in the case of systems whose characteristic families admit a 
single, connected, N — 1-dimcnsional degenerate manifold (|1.8p , it was introduced 
in [3] a decreasing potential interaction functional which is of second order w.r.t. 
the total variation (measuring the potential interaction between any pair of waves 
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as proportional to the product of their strengths, no matter if they belong to the 
same family or not). 

In the present paper, in connection with a fixed threshold parameter Sq > 0, 
we define a Glimm type functional Q = Qq + c Q, for a suitable constant c > 0, 
which is the sum of a quadratic term Qq and of the cubic interaction potential Q 
defined in [6] . Here, in presence of interactions between waves of the same families 
and strength smaller than Sq, Qq behaves as the interaction functional introduced 
in [3], while the decrease of Q controls the possible increase of Qq at interactions 
involving waves of the same family and strength larger than Sq. Employing this 
functional we can produce a simplified wave partition pattern whose errors are 
controlled by the total decrease of the Glimm functional in the time interval taken 
in consideration, and thus yield the error estimate (|1.13p . 

Note added. During the completion of the present paper, we have had knowl- 
edge of a contemporary different proof of the same convergence rate (jl.l3p provided 
by J. Hua, Z. Jiang and T. Yang [H], for Glimm approximations of a system (|1.14p 
satisfying the assumption (H). Their proof is obtained by using an adapted form 
of the functional introduced in [551 123] > that takes care of the errors in the wave- 
speeds at interactions between waves of the same family. 

2 Preliminaries 

Let j4 be a smooth matrix- valued map defined on a domain f2 C M^, with values 
in the set of N x N matrices. Assume that each A{u) is strictly hyperbolic and 
denote by {Ai(u), . . . , Xn{u)} C [0, 1] its eigenvalues. Since we will consider only 
solutions with small total variation that take values in a neighborhood of a compact 
set K G ^j^it is not restrictive to assume that O is bounded and that there exist 
constants Aq < • • • < Xn such that 

Xk-i < \k{u) <\k, Vu, k^l,...,N. (2.1) 

One can choose bases of right and left eigenvectors rk{u), lk{u), [k — 1, . . . , N), 
associated to Afe(u), normalized so that 

|rfe(^.)| = l, {h{u),rk{u))^\^^ t k^h' ^^-^^ 

By the strict hyperbolicity of the system, in the conservative case (II. ip (where 
A{u) = DF(u)), for every fixed uq e and for each k G {1, . . . , A^j-th characteris- 
tic family one can construct in a neighborhood of uq a one-parameter smooth curve 
•Sfeiuo] passing through uq (called the k-th Hugoniot curve issuing from uq), whose 
points u G >5'fc[uo] satisfy the Rankine Hugoniot equation F{u)~F{uo) — cr(u — mq) 
for some scalar a = (Tk[uo, u]. The curve S'a;[uo] is tangent at uq to the right eigen- 
vector rfe(uo) of A{uo) associated to Xk{uo), and we say that (u^, u^) is a shock 
discontinuity of the k-th family with speed (Tfc[u^, u^] if £ Sk[u''"]. 
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We describe here the general method introduced in [71 [S] to construct the 
self-similar solution of a Riemann problem for a strictly hyperbolic quasilinear 
system (jl.l4l) . As customary, the basic step consists in constructing the elementary 
curve of the k-th family [k — 1, . . . , iV) for every given left state , which is 
a one parameter curve of right states s i-^ Tk[u^]{s) with the property that the 
Riemann problem having initial data (w^, u^), = Tk[u^]{s) , admits a vanishing 
viscosity solution consisting only of elementary waves of the fc-th characteristic 
family. Such a curve is constructed by looking at the fixed point of a suitable 
contractive transformation associated to a smooth manifold of viscous traveling 
profiles for the parabolic system with unit viscosity (1.16)i. 

Given a fixed state wq € 51, and an index k € {!,..■ ,N}, in connection with 
the N + 2-dimensional smooth manifold of bounded traveling profiles of (1.16)i 
with speed close to Afc(ito), one can define on a neighborhood of {uq, 0, Xk{uo)) G 
M.^ X R X R suitable smooth vector functions (u, Vk,o-) i-^ rk{u, Vk, a) that satisfy 

rfc(Mo,0,(T) = rfc(uo) V cr , (2.3) 

and are normalized so that 

{lk{uo),rk{u,Vk,a)) = 1 y u , Vk , a . (2.4) 

The vector valued map rk{u, Vk, cr) is called the fc-th generalized eigenvector of the 
matrix A{u), associated to the generalized eigenvalue 

Afc(u,Wfe,cr) = {ikiuo), A(u)rk{u,Vk,(j)) , 

that satisfies the identity 

Afe (uo, Ufc, cr) = Afc(uo) yvk,o-. (2.5) 

Next, given a left state in a neighborhood of uq and < s << 1, consider the 
integral system 

Jo 

' Vkir) = Fk{T; u,Vk,a) - conv[o^s]Fk{T; u,Vk,<j) , 0<t<s, (2.6) 

^'^i'^) = ^conv[o,s]-Ffc(r; u,Wfe,cr) , 

where r i— > /^(r) ^ Fk{T; u,Vk,(j) is the "reduced flux function" associated to 
(1.16) defined, by 

Mr) = r M^iO, VkiO, ^(0) , (2.7) 
Jo 
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and we let conv[o,s] /a:(''') denote the lower convex envelope of fk on [0, s], i.e. 

conv[o,,] 7fc(T) = mf{e fk{y) + (1 - 0) fk{z) : 

9 e [0, 1] , J/, z G [0, s], T = 0y + {1- 0)z} . (2.8) 

Relying on (|2.3p . (|2.5p it is shown in [HIS] that, for s sufficiently small, the transfor- 
mation defined by the right-hand side of (12. 6p maps a domain of continuous curves 
r H- » (m(t), i;fc(T), cr(T)) into itself, and is a contraction w.r.t. a suitable weighted 
norm. Hence, for every in a neighborhood i^o of uq, the transformation defined 
by (|2.6p admits a unique fixed point 

r i-> (w(t; w-^, s), i;fc(r; if^, s), (t(t; u"^, s)) t G [0, s] , 

which provides a Lipschitz continuous solution to the integral system (|2.6p . The 
elementary curve of right states of the fc-th family issuing from is then defined 
as the terminal value at r = s of the u-component of the solution to the integral 
system (|2.6p . i.e. by setting 

rfe[u^](s) = u(s;u^,s). (2.9) 

Sometimes, the value (|2.9p of the elementary curve issuing from will be equiv- 
alently written Tk{s)[u^]. In the following it will be convenient to adopt the 
notations 

o-fc[u-^](s,T) = cr(T; u^, s) 

V r e [0, s] , 

Fk[u^]{s,T) = Ffe(T; u{ - ; u^,s), Vk{ - ; w^,s), cr( • ; u^,s)) 

(2.10) 

for the cr-component of the solution to (|2.6p , and for the reduced flux evaluated in 
connection with such a solution. 

For negative values s < 0, \s\ << 1, one replaces in (|2.6p the lower convex en- 
velope of Fk on the interval [0, s] with its upper concave envelope on [s, 0] (defined 
in analogous way as (|2.8p ). and then constructs the curve Tk[u^] and the map 
exactly in the same way as above looking at the solution of the integral 
system ()2.6p on the interval [s,0]. The elementary curve Tk[u^] and the wave- 
speed map ak[u^] constructed in this way enjoy the properties stated in in the 
following theorem, where we let C/([a, b]) denote the set of continuous and increas- 
ing scalar functions defined on an interval [a, 6], and we set C/([a, 6]) == C£)([6, a]) 
in the case a > b, letting Cd([5, a]) denote the set of continuous and decreasing 
scalar functions defined on [6, a] . 

Theorem 2.1 ([7l[5]) Let A he a smooth, matrix valued map defined from a do- 
main C into M^^^(IR), and assume that the matrices A{u) are strictly 
hyperbolic. Then, for every u € Vl, there exist N Lipschitz continuous curves 
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Tk[u]{s) G n satisfying lim -j^Tk[u]{s) — rk{u), together with N continuous 



functions s — )■ (Tfc[u](s, •) G C/([0, s]) (fc = 1, . . . , N), defined on a neighborhood of 
zero, so that the following holds. Whenever G fl, — Tfc[u^](s), for some s, 
letting! = {r G [0, s] : crfe[u^](s, r) 7^ (Jk[u'^]{s, t') for all r' 7^ r}, the piecewise 
continuous function 



if x/t < (7k[u^]{s, 0) , 
u{t,x) ^ {Tk[u^]{T) if x/t^ak[u^]{s,T) for some t el , (2.11) 

if x/t > ak[u^]{s, s), 

provides the unique vanishing viscosity solution {determined by the parabolic ap- 
proximation (1.16)) of the Riemann problem ([7777]), 

Remark 2.2 If the system Jj.J^I ) is in conservation form, i.e. in the case where 
A{u) — DF{u) for some smooth flux function F, and if the characteristic fields 
satisfy the assumption (H), the general solution of the Riemann problem pro- 
vided by h2.11\) is a composed wave of the k-th family made of a finite number 
of contact- discontinuities (which satisfy the Liu admissibility condition of Def- 
inition adjacent to rarefaction waves. Namely, the regions where the Vk- 
component of the solution to \2. b]) vanishes correspond to rarefaction waves if 
the a-component is strictly increasing and to contact discontinuities if the a- 
component is constant, while the regions where the Vk-component of the solution 
to \2. 6]) is different from zero correspond to contact discontinuities or to com- 
pressive shocks. In particular, whenever the solution of a Riemann problem with 
initial data u^, — Tfc[u"^](s) contains a Liu admissible shock joining, say, two 
states Tk[u^]{s'),Tk[u^]{s"), s',s" G [0,s], one has ak[u^]{s, s') = ak[u^]{s, t) 
for all T G [s',s"], and ak\u^]{s, s') provides the shock speed of the discontinuity 
{Tk[u^]{s'), Tk[u^]{s")') . Clearly, in a non conservative setting, "admissibility" 
for a jump means precisely that the jump corresponds to a traveling profile for the 
parabolic approximation with identity viscosity matrix (1.16)i. 

Once we have constructed the elementary curves Tk for each /c-th characteristic 
family, the vanishing viscosity solution of a general Riemann problem for (jl.ip is 
then obtained by a standard procedure observing that the composite mapping 

[st,...,SN)^TN{sN)o---oTi{.si)[u^]^u^ , (2.12) 

is one-to-one from a neighborhood of the origin in onto a neighborhood of . 
This is a consequence of the fact that the curves Tk [u] are tangent to (u) at zero 
(cfr. Theorem l2.ip . and then follows by applying a version of the implicit function 
theorem valid for Lipschitz continuous maps. Therefore, we can uniquely deter- 
mine intermediate states = loq, loi, . . . , oj^ = 1 and wave sizes si, . . . , sn, 
such that there holds 

uok=Tk[LOk-i]{sk) k = l,...,N, (2.13) 
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provided that the left and right states u^,u^ are sufficiently close to each other. 
Each Riemann problem with initial data 

- ( \ J^fe-i if 2; < 0, 

Uk{x) = < (2.14) 
I cjfe if a; > 0, 

admits a vanishing viscosity solution of total size Sk, containing a sequence of 
rarefactions and Liu admissible discontinuities of the fc-th family. Then, because 
of the uniform strict hyperbolicity assumption p.ip . the general solution of the 
Riemann Problem with initial data (u^, m^) is obtained by piecing together the 
vanishing viscosity solutions of the elementary Riemann problems (|l.ip (|2.14p . 
Throughout the paper, with a slight abuse of notation, we shall often call s a wave 
of (total) size s, and, if — Tk[u^]{s), we will say that (u^, u^) is a wave of size 
s of the k-th characteristic family. 

A fundamental ingredient in order to get a convergence rate for the Glimm 
scheme is the wave tracing procedure, which was first introduced by T.P. Liu in his 
celebrated paper [19j for systems with genuinely nonlinear or linearly degenerate 
fields, and lately extended to systems fulfilling assumption (H) [521 [13|. In this 
spirit, we introduce the following notion of partition of a /c-wave {u^ , u^), defined 
in terms of the elementary curves Tk at 



Definition 2.3 Given a pair of states u^,u^, with — Tk[u'^]{s) for some 
s > 0, we say that a set {y^, . . . ,y^^ is a partition of the k-th wave {u^,u^) at 
time is, if the fallowings holds. 

1. There exist scalars s^ > 0, h = 1, . . . ,1, such that, setting = ■^^^ 

= Tk[u^]{T^), there holds 

= w'' - w'^-i V h . 
The quantity s^ is called the size of the elementary wave y^ . 

2. Letting a = CTfc[u^](s, •) he the map in (j2.10p . there holds 

<j[sh) - <y{sh-i) < e V /i. 

Moreover, we require that e,+i ^]ct(t''-1), cr(T'')[, for all h (so to avoid 
further partitions of y'^ at t = {i + l)e). 

The definition is entirely similar in the case u — Tk[u^]{s), with s < 0. In 
connection with a partition {y^, . . . , y^} of {u^ , u^), we define the corresponding 
speed of the elementary wave y^ as 

^k^^f <^ir)dr yh. (2.15) 
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3 The case of a single linearly degenerate mani- 
fold 



In this section we will establish the basic estimates on the change in size and speeds 
of the elementary waves of an approximate solution provided by the Glimm scheme, 
under the following simplified assumption for the hyperbolic system (or for 

the quasilinear sytem ()1.14p ). 

(HI) For each /c e {1, . . . , 7V}-th characteristic family the linearly degenerate man- 
ifold A4k But p.8p is either empty (GNL), or it is the whole space (LD), or 
it consists of a single smooth, A^— 1-dimensional, connected, manifold and 
there holds (NGNL). 

The general solution of a Riemann problem for a sysytem satisfying the assump- 
tion (HI) consists of rarefaction waves, compressive shock and composed waves 
made of a single one-side contact discontinuity adjacent to a rarefaction wave. For 
such systems, we may consider the same type of quadratic interaction potential 
introduced in [S] for approximate solutions constructed by a front tracking algo- 
rithm, which in the case of solutions u'^ generated by a Glimm scheme can be 
defined by setting 



Co 



E + E 



(3.1) 



\SaSi3\ 



where cq > 2 is a suitable large constant to be defined later, Sa denotes the 
size of a wave of the ka-th family of u'^{t) located at Xa{t), while s'^, are, 
respectively, the (possibly zero) rarefaction and shock components of a wave Sa- 
The presence of the factor 2 in the first two summands guarantees the invariance 
of Qi when two portions of rarefaction fans of the same family, emanating from 
two consecutive mesh-points, are joined together for the effect of sampling, since 
otherwise the quantity Qi would increase for the presence of the square of the 
rarefaction components. As customary, we shall define the total strength of waves 
in u'^{t) as 



V{t) 



(3.2) 



To fix the ideas, assume that the second derivative of Afe in ()1.9p is negative, 
i.e. that 

V(VAfcTfe)(M)-rfe(ii) <0 yueMk- (3.3) 
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In order to control the nonlinear coupling of waves of the same family and with 
the same sign of two Riemann solutions for sysytems satisfying the assumption 
(HI), as in [3] we introduce the following definition of quantity of interaction. 

Definition 3.1 Consider two nearby waves of sizes s',s" with the same sign and 
belonging to the the same k-th characteristic family, with s' located at the left of 
s" . Let u' ,u" he the left state of s',s" , respectively, and assume that there exist 
waves s[, k < i < N , of the i-th family, s", 1 < j < ^, of the j-th family, so 
that u" = ( O'^ZI Tj{sJ)) o ( Oilfe T,{s'i))[u']. Then, we define the quantity of 
interaction between s' and s" as 

h{s',s")^{\{s' + s"Y-s"-\ + \s'^\)\s"^, (3.4) 

where s' + s" must be interpreted as the size of a k-wave having left state = 
Oj=i?j"(s" )[■"']; while s^ , s^ denote, respectively, the (possibly zero) rarefaction 
and shock components of a wave s. 

Remark 3.2 In the case where s' , s" are both rarefactions the quantity of inter- 
action Ii in p.4p vanishes, while /i(s',s") — \s's"\ whenever s's" are both shock 
waves. 



By standard arguments (e.g, see [HI Section 9.6, Section 13.4]) one can ob- 
tain as in [3] the basic estimates on the change in values of the total strength of 
waves V{t) and of the interaction potential Qi{t), across the grid-times ie, for an 
approximate solution u"^ constructed by the Glimm scheme. Namely, defining for 
every pair of waves of the same family s' , s" the amount of cancellation C(s', s") 
as 

fmin||s'|,|s"|| if s's"<0, 
C{s',s") = \ ^' ' ' (3.5) 
[ otherwise, 

the following generalization of [3', Lemma 2.1, Lemma 5.1] hold. 

Lemma 3.3 Under the assumption (HI), let s'l, . . . , s'j^ and s",...,s'^ be, re- 
spectively, the sizes of the waves in the solution of two adjacent Riemann problems 
{u^,u^) and {u^ ,u^), and s'( belonging to the i-th characteristic family. Call 
si, . . . , sn the sizes of the waves in the solution of the Riemann problem (u^ , u^), 
Si belonging to the i-th characteristic family. Then, there holds 



N 



|sfc-Sfc-Sfcl 0(1)- 



l<i,j<N i=l,...,N 4=1,..., AT 



(3.6) 



Moreover, for any k G {1, . . . ,N}-th NGNL characteristic family, the following 
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estimates on the rarefaction components of the outgoing waves hold. 



4 + 47 



Oil) 



E I 

. l<ij<A' 
i>j 



if 44 >o 



(3.7) 



\sl-{st + s'n\^0{l) 



E 

L l<i,3<N 
i>j 



41} 



min{l4hPfc 
„' „" 



if 44 <0: 



(3.8) 



where, in Jg. 7\ ) s'^ + s'^ represents the size of a k-wave having left state 

(0,ti'T,(4))o(oti'r,(sO)[«^]. 



Lemma 3.4 In the same setting of Lemma \3.3\ provided that the total strength 
of waves is sufficiently small, there exists some constant cq > (in p.ip ) so that 
there holds 



l<i<N 



E i^^4i + E ^i(^-4) 

l<i,j<N l<i<N 
i>3 s'-s'-'>Q 



AQi < -- 



E i^^4i + E M^4i + E ^(^-4) 

l<i,j<N l<i<N l<i<N 

i>j s'-s'-'<0 s'-s'/>0 



(3.9) 



(3.10) 



Here, as customary, we use the notations AV = — V , AQi = — , 
where V~,Qi and V~^,Qf denote, respectively, the values ofV,Qi related to the 



incoming waves Si, 



,s'^, and to the outgoing waves si 



, SAT- 



Relying on Lemma 13.41 one deduces that there exists some constant Ci > 0, 
independent of e, so that if V{t), Qi{t) denote the total strength of waves and the 
interaction potential of an approximate solution u'^{t) constructed by the Glimm 
scheme, the functional 



t^Ti{t)^V{t) + CiQi{t) (3.11) 

is non increasing at any time, provided that the total initial strength V{0) is 
sufficiently small. Moreover, for any given < m < n, the total amount of wave 
interaction and cancellation taking place in the time interval [me, ne] is bounded 
by 0(1) • (Ti(7Ti£) — Ti(n£)). Denote AT™'" ^ Ti(ri£) — Ti{me) the variation of 
Ti on [me, ne]. 
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A basic ingredient of the strategy followed in |TU] to establish a convergence 
rate of the Glimm scheme is the wave tracing algorithm introduced in (22] for 
GNL or LD systems, and then extended in [53] to NGNL systems, which consists 
in partioning the outgoing waves issuing from every mesh point {ie, je) in two type 
of waves: primary waves (i.e. waves that can be traced back from the time t — ie to 
a previous time t = me < ie), and secondary waves (i.e. waves that are generated 
by interactions occurring in the time interval ]me,ie], or that are canceled before 
a later time t = ne > ie). The total strength of secondary waves produced in 
a given time interval [me, ne] is bounded by the total amount of interaction and 
cancellation occurring within [me,ne]. 

The key step of this procedure is to show that the variation of a Glimm func- 
tional provides a bound for the change in strength and for the product of strength 
times the variation in speeds of the primary waves. The main novelty of the analy- 
sis performed here consists in implementing a wave tracing algorithm for a NGNL 
system satisfying the assumption (HI) in which such bounds are obtained relying 
on a Glimm functional with a quadratic potential interaction, differently from the 
Glimm functional with a cubic potential interaction used in [23' . Namely, recalling 
the Definition 12.31 of a wave partition, we have the following result. 

Proposition 3.5 Under the assumption (HI), given a Glimm approximate solu- 
tion and any fixed < m < n, there exists a partition of elementary wave sizes 
and speeds {y'i!:{i,j), A['(i, j)}, k = 1, . . . , N , i — m,m + 1, . . . ,n, j £ Z, so that 
the following hold. 

1. For every i,j,k, {2/fc (*; i)}o<?i<£^(i j) pO'i"t'ition of the wave of the k-th 
family issuing from {ie,je), and {^ki''-' j)} Q<:h<ik{i j) corresponding 
.speeds, according with Definition \2.S\ 

2. For every i,j,k, {yl{i,3), ^fc(«7 j)}o<ft<£^(jj) « disjoint union of the two 
sets 

[vl^hj), A^(i,j)} , {yk(.h3), Kih])] , 

with the following properties: 
(a) 

5^ ||y;(*,j)|| =0(l).|ATr"| Vm<i<n; (3.12) 

j,k,h 

(b) for every fixed i,k,h, there is a one-to-one correspondence between 
{yl{m,j),\l{m,j) ; j G Z} and {yUhj)Xihj) ; J G Z}- 

{y^(m,j),A^(m,j)} ^ {yUhi{^,J,k,h))Xkihei^,J,k,h))} (3.13) 
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such that the sizes and the speeds of the corresponding waves 
satisfy 



y I max 

j,k,h 



E 



- sU^Ji^,J,k,h))\^ = Oil) ■ \ATT'''\ , (3.14) 
f|s^(m,j)|- max |A^(m,j) - A^(^, = • |ATr"| ■ 

\ m<i<n / 



(3.15) 



Proof. The desired partition for an approximate solution will be constructed 
proceeding by induction on the time steps is, m < i < n. Assuming that a 
partition of elementary waves fulfilling properties 1-2 is given for all times me < 
t < ie, we wish to produce a partition of the outgoing waves generated by the 
interactions occurring at t = ie, so to preserve the properties 1-2. Observe first 
that the existence of such a partition is already guaranteed by the analysis in 
[23j if all interactions take place between waves of different family or of the same 
family with opposite sign, since for systems satisfying the assumption (HI) the 
change in strength and the product of strength times the variation in speeds of the 
primary waves is controlled by the variation of a Glimm functional with quadratic 
interaction potential as the part in brackets of (|3.ip . 

Therefore, it will be sufficient to consider an interaction between two waves 
issuing from two consecutive mesh points ((i — l)e, (j — l)e) and ((z — 1)£, je), say 
sjj,, s'l, belonging to a k-th NGNL characteristic family, and having the same sign. 
For the sake of simplicity, assume that sj^,, s'^' > 0. Let Sp [p = 1, . . . ,N) be the 
outgoing wave of the p-th family issuing from {ie,je), and let 

{y'k, ^'k}o<h<e' ' {yk''^ '^k'''}o<h<e" ' {3-16) 

be the partitions of sJc and s'^ enjoing the properties 1-2 on [me, {i — l)e], with 
sizes 

0<h<£' ' Q<h<e" ■ 

(3.17) 

For every p ^ fc-th wave Sp, we may choose a partition {yp}Q<h<ep sls in Def- 
inition 12. 3[ with corresponding speeds {^p}o<h<ip- Then, if we label all the 

subwaves as secondary waves j/p, the bound p.l2p (for i,j,p) is certainly 
satisfied thanks to the interaction estimates (|3.6p . Instead, for the A;-th wave 
Sfc, possibly considering a refinement of the partition of s'^ (or of s'f.) we may 
assume that either s'j, -I- s'^^ < Sk, or s'^ < Sk (in the case s'j^ > Sk), and let 

/ = max{/i < £' : J2g=i 4 < Sk}, t = max{/i < i" : 4 + Y!1^^ 4" < Sk}- 
Then, we define a partition of s^ by means of its sizes, setting 



4'' ii h^i,...,i , 

(3.18) 

sl''-*^' if ^' = £' and £' + !,...,£' + /' 
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(possibly refining the partitions p.l6p so to satisfy property 2 of Definition I2.3p . 
and choosing a partition oi Sk — {s'f.+s'^) as in Dcfinition l2.3l in the case Sk > Sfe+s^'- 
The subwaves in p.lSp inherit the same classification in primary and secondary 
waves of the corresponding subwaves s'^ or s^."*"^ , while all the possible subwaves 
of Sk — (sfc + s'l) are labelled as secondary waves. Clearly, the bound (|3.12|) is 
again satisfied because of the interaction estimates (|3.6I) . while the ono-to-one 
correspondence at (|3.13p and the bound p.l4p are verified by construction and 
by the inductive assumption. Hence, in order to conclude the proof, it remains to 
establish only the estimate p.lSp . 

By the assumption (HI), and because the incoming waves sj^, s'^ have the same 
sign, at most one of them can possibly be a composed wave, say s'^, while s^' 
will be a shock. Denote as {s'/^Y, (s^)'* the rarefaction and shockcomponent of 
sjj., respectively. For sake of simplicity, assume that Sk > s'j,, i.e. that £ — £'. 
The outgoing wave Sk is either a shock or a composed wave. In the first case 
its Rankine-Hugoniot speed Afe coincides with the speeds of all subwaves 
defined according with Definition l2.31 since for a shock wave the integrand function 
in (|2.15p results to be a constant (cfr. Remark 12. 2p . Hence, letting (A'^)*, A'^' 
denote the speeds of the shock component of sjj and of s'l, respectively, by a direct 
computation one finds 



lA^-A'/'l < \Xk-{X'kr\=0{l)-4 \fh^ 

|A^ - Af I = |Afc - A',:| = 0{1) -4 V h = £' + !,...,£' +1" 

In turn, p.l9p implies 



(3.19) 



h=l h=t' + l 

which, relying on the inductive assumption, yields (|3.15p since in this case, by the 
estimate p.7p . and because the rarefaction component of Sk is zero, there holds 

Next, assume that the outgoing wave Sk is made of a rarefaction component 
{skY < (sfc)'' and of a shock component {skY ■ Then, possibly considering a 
refinement of the partition of Sk (and hence of the partition of sj,), there will be 

some index t < £' so that Yl\i=i ^k — i^kY ^.nd Yl\i=i^+i ^k = {^kY ■ Notice that 
{skY can be seen as a shock wave generated by an interaction between a composed 
wave with rarefaction component of size {4Y ~ i^kY ^^d shock component {s'i.Y, 
and of the shock wave sj!, for which we can apply the above estimates on the 
variation of wave speeds. Hence, the wave speeds Aj^ of sj! defined as in (|2.15p 
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satisfy 



fo v/i = i,...,r, 

- Af I = . [(4)'- - + (4)1 yh = e' + i,...,e' + I 



This implies 



(3.21) 



e' e'+e" 

5^4'|A^-A11+ E 4'1A^-Af |=0(l).[((4)'--(.fen + (4r]4', (3.22) 

which in turn, relying on the inductive assumption, yields again p.lSp since in 
this case, by the estimate p.7p and because s'l = {s'D^, one has [(sj.)'" — {skY + 
(■^fc)"] ■5fc = '^(l) ' -^i('*fc' ^k)- This completes the proof of the proposition. □ 

Proposition l3.5l provides for NGNL systems satisfying the assumption (HI) the 
same type of result that was established in (TUl Proposition 2] for systems with 
GNL or LD characteristic families. In order to obtain the desired convergence 
rate p.l3p one can now simply repeat the proofs of [101 Propositions 3-4] and 
of the final estimates in [101 § 6], which all rely only on the conclusion of [lOl 
Proposition 2] and thus remain valid within our more general framework of NGNL 
systems. We will give a brief description of them in Section [51 

Remark 3.6 The conclusion of Theorem \1.2\ established so far for smooth sys- 
tems satisfying the assumption (HI), remains valid if we assume that the flux 
function F is C^'^ and that, for each k-th characteristic family not fulfillying (HI), 
the linearly degenerate manifold M.k in (II. 8p is a C^'^ N -dimensional, connected 
manifold, F is on fl \ Aik, the vector field rk is transversal to the boundary of 
Mk, dnd letting d~^A4k, d~A4k denote the connected components of the boundary 
of Mk where points towards il,\Mk and Mk, respectively, there holds 

V+(VAfc-rfc)(u)-rfc(M) <0 Wed+Mk, 

(V±(VAfc • rk){u) ■ rk{u) ^ ^ lim^ ^"+''"'° (^"))-^^'° ^''^ denoting the one-side 

second derivatives of Xk)- Indeed, the only difference in the structure of the ele- 
mentary waves of a NGNL k-th family satisfying such assumptions instead of (HI ) 
comes from the possible presence of two-sided contact discontinuities. In fact, un- 
der the above assumptions, the general solution of a Riemann problem of the k-th 
family will be either a rarefaction wave, or a shock wave (which can be either a 
compressive shock or a contact discontinuity), or a composed wave made of a rar- 
efaction wave adjacent to one (one-sided or two-sided) contact discontinuity or 
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several (two-sided) contact discontinuities. Then, we may consider the interaction 
potential Qi in (|3.ip . where the shock component of a composed wave Sa containing 
several contact discontinuities j^, ■ • ■ , / = X]p=i ^a,p! and for every such 
wave we add the term 2 ^ p,? |Sq „ „|. One can easily verify that employing this 

definition of Qi and the same definition of quantity of interaction Ii in p.4p . the 
estimates stated in Lemma \3.3\ continue to hold, provided that 

inf |s > : Rk[u]is) e d+Mk , u e d'Adk^ > (3.24) 

(Rk[u\{s) denoting the integral curve of rk), which is certainly true up to a pos- 
sible slight restriction of the domain f2. Relying on Lemma \3.3l one then deduces 
Lemma \3.4\ and thus can establish the key Provosition l 3. 5\ with the same arguments 
as above. 

4 A new interaction potential 

We turn now our attention to an approximate solution u"^ constructed by the 
Glimm scheme for an hyperboUc system (|l.ip (or for the quasihnear system (jl.l4p ) 
that satisfy the assumption (H) stated in the Introduction. We recall that for such 
systems the general solution of a Riemann problem contains composed waves made 
of several contact discontinuities adjacent to rarefaction waves (instead of just 
a single contact discontinuity adjacent to a rarefaction wave as for the systems 
treated in § We will say that a wave s of this type, belonging to the k-th 
characteristic family, crosses all connected components of A4k that are transversal 
to the fc-th elementary curve Tk issuing from the left state of s and terminating 
on the right state of s. Notice that, for each fc-th NGNL family, and for every 
connected component MkM of Mk, the first derivative VAfc • rk has opposite signs 
on the connected components of fl\A4k.h adjacent to Mk,h, and as a consequence 
the second derivative V(VAfc - rk) - rk has opposite signs on any pair of consecutive 
components A4k,h,Mk.h-^'l■ Thus, by continuity we may assume that there exists 
some constant (5o > so that 

V{VXk-rk)iu)-rk{u)^0 y u s.t. d^u, Mk) < 6Sq , (4.1) 

where d{u, A4k) ^'oiiveMk I"" ~ ''^1 denotes the distance of a state u from Mk- 

Remark 4.1 Conditioner^ implies that every wave s of a k-th NGNL family 
with strength \s\ < 3Sq crosses at most one connected component of Mk- Moreover, 
if an interaction takes place between two waves of the k-th characteristic family 
with strength < Sq, then, by the interaction estimates in |3 Theorem 3.7], the 
outgoing wave of the k-th family crosses as well at most one connected component 
ofMk. 

By Remark |4. 11 as far as the waves of the NGNL families involved in an interac- 
tion have all strength smaller than Sq, we can establish the same kind of estimates 
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of Proposition 13.51 employing the quadratic interaction potential in (|3.ip even for 
systems satisfying the more general assumption (H). On the other hand, observe 
that if we consider an interaction between two shock waves of a k-th NGNL fam- 
ily, say s', s", with speeds A', A", respectively, and we assume that s' , s" have the 
same sign, then, letting A denote the shock speed of the outgoing wave of the fc-th 
family, by the interaction estimates in Theorem 3.7] there holds 



|s's 


"1 


|A'- 


X"\ 




s' 


+ s" 





[sAX] = |.'||A - A'l + |."||A - A"| = 0{1) ■ ' I ' „| ' . (4.2) 



Notice that |s's"| | A' — A"| has precisely the same order of the quantity of which 
it decreases the interaction potential Q introduced in [6] whenever interactions of 
this type take place. Therefore, if we assume that at least one of the incoming 
waves of the k-th family has strength > So, we deduce from (|4.2p that [sAX] = 
0(1) • |AQ|/(5o. Hence, for such interactions one may derive the same kind of 
estimates on the products of the wave strengths times the variation of the wave 
speeds of Proposition 13.51 employing the cubic interaction potential Q defined in 

In view of the above observations, we shall introduce now a functional Q that 
is the sum of a quadratic and of a cubic interaction potential. The latter is the 
interaction potential for waves of the same family and with the same sign defined 
in [6], valid for general strictly hyperbolic systems (|1.14p . which takes the form 



,, Jo Jo 



(4.3) 



The summation here extends to all pair of waves Sa, sp of the ka € {1, . . . , N} 
family with the same sign (including Sq — sp), of the approximate solution u'^(t), 
and ctq (Tfe^ [cJq,](sq, •) denotes the map in (j2.10p . where cJq, is the left state of 
Sa- Such a functional controls the nonlinear coupling of waves of the same family 
with the same sign. 

The quadratic part Qq of the functional Q enjoys two basic properties: 

1. it decreases whenever it takes place an interactions between "small" waves 
of the same family, i.e. waves whose strength is smaller than and the 
amount of decreasing satisfies the same type of estimate p.lOp obtained for 
systems with a single linearly degenerate manifold; 

2. the possible increase of Qq caused by interactions involving "large" waves of 
the same family, i.e. waves of strength larger than Sq, is controlled by the 
decrease of Q. 

Thus, for general hyperbolic systems (|1.14p satisfying the assumption (H), we shall 
consider a potential interaction of the form 

Q{t)^Qg{t)+cQit), (4.4) 
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where c > 2 is a suitable constant to be specified later. 

Towards the defintion of Qq, let us first introduce some further notations. Given 
a composed wave s of a fc-th NGNL family, let {s'^}h=i,...,i be its decomposition in 
rarefaction and shock components, and write h G TZ (respectively /i G 5) if s'' is a 
rarefaction (respectively a shock) wave. Thus, letting w'^~^, w''- denote the left and 
right states of each wave s'', one has — T/j[w''~^] (s''). Next, for every given 
shock s^, h G S, we define a convex-concave sub decomposition {s'''^}p=i_...^g^ 
as follows. Assuming for the sake of simplicity that s'' > 0, let — < < 
■ ■ ■ < t'"^ — he & partition of determined by the inflection points of the 
reduced flux r ^ fj^ir) ^ Ffc[u;''-i](s'', r) in^ (ITTO)l . and set s'^^p = tP - tP-^. 
We will write p G ^ (respectively p € '~-) if is convex (respectively concave) 
on [tP 



fP], and we will call s'^'P a convex (respectively concave) component of 



s if p G ' (respectively p € Then, considering the affine map 



1 if \s\>25o, 

cp{s)^{ {\s\-So)/So if So< \s\<2So, 
if |s|<<5o, 

we define the intrinsic interaction potential of s'*, h £ S, as 



(4.5) 



qis") = vis") 



„h,p\ 



(4.6) 



where the first summand runs over all indexes p,q U p ^ q, and q{s'^) 

is understood to be zero if has zero convex component. Notice that, by defi- 
nition ()4.5p . for shocks s'* with non zero convex compoents, q(s'') can possibly be 
zero only when h = 1 or h = I, i.e. when is the first or the last component of 
s. In fact, all other shock components of s are two-sided contact discontinuities 
which necessarily must cross at least two connected components of A^fc, and hence 
their strengths are certainly larger than 2do because of ()4.ip . 

Now, defining the inner interaction potential of a composed wave s as 



hes,KeTZ 



E 

hes 



El 



(4.7) 



we can finally provide the definition of the quadratic interaction potential enjoing 
properties 1-2 by setting 



Qg{t)^2 \saSp\ + Y,Q'isc.)+c 



E + E 



SqS/3 



(4.8) 



where, as usual, Xa{t) denotes the position of the wave Sq,, and ka its characteristic 
family while c is the same constant that appearzs in ()4.4p . Here, the second 
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summation runs over all composed waves Sa present in u'^ (t) . Notice that Qq differs 
from the interaction potential Qi defined in § [3] only for the presence of the inner 
interaction potential of the composed waves that replaces the corresponding 
terms of the second and third summands in On the other hand, whenever 

|sa| < So, we clearly have Q^(sq) = kaS^I + |sqP, h e S, k £ TZ, and thus one 
recovers the same expression present in Qi. 

Remark 4.2 Consider a shock wave s with strength \s\ < Sq that crosses a con- 
nected component Mk,h of Mk- According with the above definitions s is de- 
composed in a convex and a concave component s^ , s"". Relying on Proposi- 
tions 2.1-2.2], we deduce that, choosing 5q sufficiently small, there holds 



< ci Is" 



(4.9) 



for some constant < ci < 1. Such a bound will be useful in the study of the 
variation of the intrinsic interaction potential q{s) in presence of interactions. 
Notice that the above estimate holds even in the case, instead of ()1.9|) . we assume 
that there is some even index p so that the following weaker condition is satisfied: 



k,h ■ 



(4.10) 



where Dl.^Xk{u) denotes the j-th derivative of Xk along rk, inductively defined by 
setting Dr,^Xk{u) ^ VAfc(ii) • rk{u), and D^^Xkiu) = Dl'^^Xk{u) ■ rk{u) for all 
J>1- 

Towards an analysis of the interaction potential above introduced, we first de- 
fine a quadratic quantity of interaction as in Section [3] for waves of the same family 
and with the same sign, to measure the decrease of the quadratic functional Qq in 
when waves of this type with strength < 5q are involved in an interaction. 



Definition 4.3 Consider two nearby waves of sizes s' , s" with the same sign and 
belonging to the the same k-th characteristic family, with s' located at the left of s" . 
Assume that \s'\, \s"\ < Sq and, with the same notations of Definition \3.i[ suppose 
that the state belongs to the connected component of^l\A4k lying between two 
consecutive manifolds Mk.h-i, Mk.h of Mk- Then, in the case s' , s" > 0, we 
define the quantity of interaction between s' and s" as 



/(s',s") = 



' [\is' + s"Y - s' 



f |s'1]|s"^| 

if V(VAfc •rfc)(w) •rfc(u) 



tf V(VAfc •rfc)(w) •rfe(w) 



<o, 



> 0. 



(4.11) 



' < 0. For notational con- 



An entirely similar definition is given in the case s' , s 
venience we also set I{s' , s") = for every pair of waves s' , s" of the same family 
that have opposite sign. 
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Next, following [6, Definition 3.5], we introduce a definition of quantity of 
interaction for a general strictly hyperbolic system (|1.14p . which measures the 
decrease of the cubic functional Q in (|4.3p when waves of the same family and 
with the same sign interact together. 

Definition 4.4 Consider two nearby waves of sizes s', s" with the same sign and 
belonging to the the^same k-th characteristic^family, with left states u' ,u" , respec- 
tively. Let F" = Fk[u']{s', •) and F" = Fk[u"]{s" , ■) be the reduced flux with 
starting point u' , u" , evaluated along the solution of 112. 6]) on the interval [0, s'], 
and [0, s"], respectively (cfr. def. h2.10\) ). Then, assuming that s' > 0, we say that 
the amount of interaction J(^s' , s") between s' and s" is the quantity 



Jis\ s") - / conv[o,,.]i^'(0 - conv[o,s'+s"]F'UF"{0 
Jo 



(4.12) 



F'is') + conv[o,,,,]F"(e - s') - cony[o,s'+s"]F'UF"{0 



d^, 



where F'UF" is the function defined on [0, s' + s"] as 



~ ~ (F'(s) if se[0,s'], 

F'UF"(s) = 1 ' ^ ^ ^ (4.13) 

\f'{s')+F"{s~s') if se[s\s' + s"]. 

Here, convjQ ;,]/; conC[ci^b]/ denote the lower convex envelope and the upper concave 
envelope of f on [a,b], defined as in ()2.8p . In the case where s' < 0, one replaces 
in 114.12^ the lower convex envelope with the upper concave one, and vice-versa. 
As in Definition (|4.3p . for notational convenience we also set J^{s',s") = for 
every pair of waves s', s" of the same family that have opposite sign. 

Remark 4.5 Notice that by the Lipschitz continuity of the derivative (u, s) i—f 
DrFkMis,-) of the reduced flux (j2lT))l ^cfr. J^), it follows J{s',s") = 0{l)-\s's"\. 
Moreover, by Remark \2.2\ one can easily verify that, in the conservative case, if 
s' , s" are both shocks of the k-th family that have the same sign, then the amount 
of interaction in takes the form 

i.e. it is precisely the product of the strength of the waves times the difference of 
their Rankine Hugoniot speeds. 

Relying on the results in [61 Section 3] and on Lemma 13.31 we will show now 
that the interaction potential Q defined by (|4.3p . (|4.4p . (|4.8p . is decreasing at 
every interaction, and that the variation of the total strength of waves V in an 
approximate solution u"^ is controlled by |A(5|. 
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Lemma 4.6 Under the assumption (H), in the same setting of Lemmas \3.S'\3.4\ 
there exists some constant c > (in (|4.4p . (14. 8p ). so that there holds 



Ay < 



0{1) 



l<i<N 



E 

l<i,j<N 
i>j 



l<i<N 



(4.14) 



AQ<-- 



E 

l<i,j<N 
i>j 



E 

l<i<Af 



E I 

l<i<Ar 
s's''<0 



E ^(^-^") 

l<i<Af 
Ms"|<5o/2 



(4.15) 



0(1) • V 



•E 

l<i<Ar 



Proof. A proof of the estimate (|4.14p can be found in 6 , thus wc wiU focus our 
attention on (|4.15p . For the sake of simphcity, we shaU consider only the case in 
which the two adjacent Riemann problems are solved by a single wave, say s' and 
s", s' on the left of s" . We distinguish three cases, depending on the strengths of 
s' and s" and on their characteristic families. 

1. s' and s" are waves of the k' and k" < k' characteristic families. 

To fix the ideas, let s', s" > 0. Observe that for every NGNL fc-family, since 
condition (|1.9p implies that the characteristic vector field is tranversal to 
A^fc, by construction it follows that the fc-elementary curves Tk are transver- 
sal to each manifold Mk.h- As a consequence of this property one can easily 
verify that, letting u', u" be the left states of s', s", and denoting u'^,u"^ the 
left states of the outgoing waves Sf and Sk" of the k' and k" characteristic 
families, there holds 



E" 





E E ^"'^ 


- E E 


+ 


Ken. 


Ken 












EE^'"^^ 


-EE4f 


+ 




Ken 


0{l)-\\Tk4n"] 


-Tk4n"+]\\ . 

II L°° 







(4.16) 



Here the L°° norm in the first and second equality is referred to the in- 
tervals [0, min{s', Sfc'}] and [0, min{s", Sfc//}], respectively. Then, since the 
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interaction estimates in jB] Section 3] imply 



Tk.[u'+]~Tk'[u'] 
Ty[u"+]-Tu.-[u"] 



0{1) 



(4.17) 



relying on (|4.16p - (|4.17p . we deduce the following bounds on the variation of 
the inner interaction potential at (|4.7p 



AQ^(s') = 0{\) ■ S' • s' = 0{1) ■ s's" , 
AQ\s") = 0(1) • S" • s" = 0(1) • s's" . 
Hence, using l|4.14p . (|4.18p . one obtains 

AQ < -c\s's"\ + 0(1) • \s's"\ + 0(1) • |s's"| • V- , 
from which we derive (|4.15p . choosing c > sufficiently large in 



(4.18) 



(4.19) 



2. s' and s" are both fc-waves and s's" < 0. 

By defintion of , and with the same analysis in the previous point, one 
deduces that in this case the inner interaction potential of the outgoing k- 
wave s satisfies Q^(s) < m.m{Q' {s'),Q^ {s")} + 0(1) ■ \s's"\. Hence AQ^ = 
0(1) • |s's"|, and thus we obtain the same estimate in (|4.19p . On the other 
hand, relying on [HI Proposition 4.1] we derive 



AQ<0{l)-V- -Cis'^s") 



(4.20) 



which, together with (|4.19|) . yields (|4.15p . choosing c > sufficiently large 
in (Ell). 



3. s' and s" are both fc-waves and s's" > 0. 

To fix the ideas, let s' , s" > 0, and call s the outgoing fc-wave. We shall 
distinguish a number of cases, depending on the strengths of s', s" . 

(a) max{s', s"} < 6o/2. 



In this case, by definitions (|4.5p - (|4.8p one has q{s') — q{s") = = q{s' + 
s") = 0, and AQq < Aq + AQi + 0(1) • \AV\, where AQi denotes the 
variation of the interaction potential Qi in (|3.ip (related to the waves 
involved in the interaction). Hence, relying on ()4.14p . and applying 
(I3in . we deduce Aq = 0(1) • \AV\ = 0(1) • J{s', s") and 



AQ,<-ll{s',s") + 0{l)-J{s',s"). 



(4.21) 
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On the other hand, due to [H Proposition 4.1], we get 



(4.22) 



which yields (|4.15p choosing c > sufhciently large in (|4.4|) . 
(b) 5q/2 < max{s', s"} < 2So ■ 

To fix the ideas, assume that s' crosses a connected component Aik.h 
of A^fe where V(VAfc ■ rt) ■ rk < 0. Because of ()4.ip . this imphes 
that the wave s" on the right of s' must be a shock with zero convex 
component and hence q{s") — 0. For sake of simplicity, we shall treat 
only the case in which also s' is a shock and s' < Sq < s/2, the other 
cases being similar or simpler since for such values of s', s there is the 
largest possible increase of q due to the fact that, by definitions (|4.5p . 
(|4.6p . one has q{s') = f{s') = 0, ip{s) — 1. Under these assumptions, 
by definitions gH), (g!?]) we have Q^{s') = Q^{s") = 0, and letting 
s"",s"", s''~",s''^, denote the convex and concave components of s,s', 
relying on (|4.14p we deduce 



Qi{s) =q{s) ^2s-s- + {s- f 

= 2s'- {s'- + s") + + 0{l) ■ \AV\ (4.23) 

<2s'-{s' + s")+Oil)-J{s',s"). 

Moreover, observe that s' < So < s/2 imphes s' < s" + 0(1) • J{s', s"). 
Therefore, using (|4.23p . and recalling that by Remark 14.21 one has s'^ < 
s''^ , we find 



Hence, (|4.24p together with (|4.22p . that continues to hold, yields (|4.15p 
choosing c > sufficiently large in (|4.4p . 

(c) min{s', s"} < 2So < max{s', s"}. 

To fix the ideas assume that s' is a composed wave of size s' < 26o, 
crossing a connected component of Mk where V(VAfc • r^) • rj, < 0. Be- 
cause of (|4.ip , this implies that the first component s"^ of the (possible 
composed) wave s" on the right of s' must be a shock of size s"^ > 2So. 
For sake of simplicity we shall treat only the case in which also s' is a 
shock, the other cases being similar. Observe that, letting u', u" be the 
left states of s', s" , calling m'+ the left state of the outgoing fc-wave s, 
and letting s'"",s''~' denote the convex and concave components of s'. 



AQg < -2s' s" + Q\s) + 0{l) ■ \AV\ 

< -2s"{s'- ~s'-)+0{l)-J{s',s' 
<0{1)-J{s',s"). 



(4.24) 
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by the same arguments at point 1 wc find 



+ 



E ^"^ - E 



Oil) 



\Tk[u'] - Tk[u'+]\\ + 



n[u"]-Tk[u'+]is' + -)\ 



(4.25) 

where the L°° norm of the two terms in the equahty is referred to the 
intervals [0, s'] and [0, min{s", s — s'}], respectively. Then, applying the 
interaction estimates in [&, Section 3], we derive 



(4.26) 



On the other hand, observe that by definition ()4.5p the above assump- 
tions imply <f{s"-^) — 1, (p{s^) — 0{1) ■ J{s',s"), since the first com- 
ponent of s satisfies the lower bound > s"^ -I- 0(1) • J{s' ^ s") > 
25o + 0(1) • J{s',s"). Thus, relying on gH]), and because 

s' < s"^, we obtain 



AQ' < 2s'- {s'- + s"i) + (s'-)2 + 0(1) . (E + \AV\) 
< -2s"\s'- - s'-) + 2s' s"^ + 0{\) ■ J{s', s") , 

which in turn, recalling that by Remark 

AQg < -2s' s"^ + AQ^ 
< 0(1) -{J is', s") 
= 0{1)-J{s',s"). 



one has s'"" < s'" 

0{1)-\AV\ 
HAF|) 



(4.27) 
yields 

(4.28) 



Hence, (|428l) together with (|4?22l) . that continues to hold, yields ([415]) 
choosing c > sufficiently large in (|4.4p . 

(d) min{s', s"} > 2Sq. 

We shall treat only the case in which the last componenent s'' of s' and 
the first component s"^ of s" are both shocks of size > 2So. The other 
cases are simpler or reducible to one of the previous cases (a), (b), (c). 
Then, by definition (|4.5p there holds <p(s'' ) = fis"^) = 1. Moreover, 
with the same notations and with the same arguments of point (c) , we 
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have 

E = 



E E ^'"^^ - E E 



+ 



Ken 



Ken 



= 0{l) 



\Tk[u']-Tk[u'H 



\Tk[u"]-n[u'+]{s' 



0{1)-J{s\s"). 



(4.29) 



Thus, relying on (|4.14p . (|4.29p . we derive 

AQ, < AQ^ + 0(1) • |Ay| 

= 0(1)- (S+|A1/|) (4.30) 

= 0(l)-^(s',s"), 

which, together with (|4.22p . that continues to hold, yields (|4.15p choos- 
ing c > sufficiently large in ()4.4p . 

□ 

Relying on the above result one can prove that there exists C > so that, 
assuming V^(0) sufficiently small, the Glimm fmictional 

t^T{t)^V{t)+CQ{t) (4.31) 

is non increasing at any time, and at every discrete time t = is there holds 



AT(ie) < —-^[amount of cancellation at t = ne] + \ AQ{ie) 



(4.32) 



Hence, for any given < to < ti, the total amount of wave interaction and can- 
cellation taking place in the time interval [me,ne] is bounded by 0(1) • IT™'"!, 
where 

AT™'" = T{ne) - T{me) 
denotes the variation of T on [me,ne]. 



5 Wave tracing for general non genuinely nonlin- 
ear systems 

We will show now how to implement a wave tracing algorithm for a NGNL system 
satisfying the assumption (H) so that the change in strength and the product of 
strength times the variation in speeds of the primary waves be bounded by the 
variation of the Glimm functional in (|4.3ip . Namely, recalling the Definition 12.31 
of a wave partition, wc have the following result analogous to Proposition 13.51 
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Proposition 5.1 Under the assumption (H), the same conclusions of Proposi- 
tionlMhold, with AT™'" m place o/ AT^^". 



Proof. As in the proof of Proposition 13.51 in order to produce a partition for an 
approximate solution u"^ that fulfills properties 1-2, one may proceed by induction 
on the time steps is, m < i < n. Then, assuming that such a partition is given 
for all times me < t > is, our goal is to show how to define a partition of the 
outgoing waves generated by the interactions that take place at t = ie, preserving 
the properties 1-2. As observed in the proof of Proposition [331 it will be sufficient 
to focus our attention on interactions between waves of the same family and with 
the same sign, since whenever any other interaction occurs for a system satisfying 
the assumption (H), the change in strength and the product of strength times 
the variation in speeds is controlled by the variation of a Glimm functional with a 
quadratic interaction potential as the part in brackets of (14. 8p (cfr. P?l Lemma 3.2 
and Theorem 5.1]). 

Thus, consider an interaction between two waves, say s'j,, s'^, issuing from two 
consecutive mesh points {{i — l)e, (j — l)e) and ((? — l)e,je), belonging to a fc-th 
NGNL characteristic family, and having the same sign. Observe that, if |sfc|, IsJ^'j < 
So/2, then relying on the estimates (|4.14p . (|4.15p provided by Lemma WM one 
obtains the desired partition proceeding precisely as in the proof of ProDOsition l3.5l 
Hence, we shall treat only the case where max{|s5.|, |sj,'|} > So/2. For sake of 
simplicity, we assume that sj,, s'/. > and that the outgoing fc-wave Sk is a shock, 
the other cases being entirely similar. Let 

{y'ki ^'k}o<h<e' ' {y'k^' '^k^}o<h<e" ' (^■^) 

be the partitions of s'j. and s'/. enjoing the properties 1-2 (on the interval [me, (i — 
l)e]), with sizes 

a<h<i' ' a<h<i" ■ 

(5.2) 

Then, define a partition {?/p}o<h<£p of the outgoing wave Sp {p = 1,...,N) of 
the p-th family issuing from (ie,je) (with corresponding speeds {^p}o<h<i.j,) as in 
Proposition 13. 51 In particular, a partition of Sk is defined by means of its sizes as 

if h^i....,e', 

(5.3) 

if / = and e' + !,...,£' + i" 




(with the same notations of the proof of Proposition [ 

Clearly, such partitions continue to satisfy the bounds (|3.12|) . (|3.14p and the 
one-to-one correspondence at (|3.13|) . thanks to the estimate ()4.14|) . and because of 
the inductive assumption. Therefore, in order to conclude the proof, it remains to 
establish only the estimate p.lSp on the wave speeds. To this end, notice that the 
Rankine-Hugoniot speed of the outgoing fc-wave Sk coincides with the speeds 

of all subwaves defined according with Definition [531 since for a shock wave 
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the integrand function cr(-) in (|2.15p results to be a constant (cfr. Remark I2.2p . 
Moreover, by the choice of the speeds of a partition at (|2.15p . one has 



a"{T) dr, 



(5.4) 



where t^'' 



Eh- 'P -.Ilh ^ Y^'i- jip 
p=l T 'k ~ Z^p=l 



s^^ , and 



a'(.) = afeK,](4,.), a"(.) = a,K](4,.): 



denote the map in (|2.10p defining the speed of the rarefaction and shock compo- 
nents of s'j, and s'^', respectively (a;^,u;^' being the left states of s'^jsj!). 

Then, relying on the interaction estimates in [BJ Section 3], with the same type 
of arguments used in the proof of [6, Lemma 3.9] one obtains the following estimate 
on the wave speeds, similar to the one provided by 23, Theorem 3.1] : 



Afc • (4 + 40 = 



a(T) dr + 0(l).^(4,4) 



g\t) dr + / a"(r) dr + 0{\) ■ J{sl 4) 



(5.5) 



= E ''^^'t + E + ^(1) • '^(^'^' ^^0 ■ 

h=l h=l 

Thus, since by the monotonicity property of cr'(-) and cr"(-), we have 

- 0{1) ■ J{s'„ 4) < Afe < 4^ + 0{1) ■ J{4,4) V h 
using (|5.5p we derive 

lA;,''-Afcl = A',''-Afe + 0(l).^(4,4) 

1 



0(1) 



p=i,...,e' p=i 

^(4,4) 



which, in turn, yields 



h=l 



4 + 4 



EE^M''(A;f-Ar)+ 

h=i p=i 



EE«'(^"-^r) 



h=i p=i 



+ o(i).^(4,4)-- 



(5.6) 
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Notice that the terms of the first double sum on the right hand side of (|5.6p are 
antisimmctric in {h,p), and hence the first summand vanishes. Moreover, recalling 
we have 



h=i p=i 



p=l h=l'^^k h=l p=l-^'^k (57-) 



a'{T)dT-s'J a"{v)dii 



k 

^0 



Jo 



[cr'(r) — cr"{r])\ drj dr . 



On the other hand, observe that the term in (|4.3p corresponding to the outgoing 
shock wave Sk vanishes (being the map a constant), and hence one clearly has 

f ' f ' \a'{T) ~ a"{r))\ di^dr = • \AQ{ie)\ . (5.8) 

JQ Jo 

Thus, since the assumption maxjsj,, s';,'} > So/2 imphes sJc + *fc > ^o/2, from 
(ESI)- 153 it follows 



(5.9) 



E4'|A1?'-Afc|<|- f" f " \a\T)-a"{v)\drjdr + Oil)-J{.sls'l) 



< 0(1) • \AQiie)\ 



An entirely similar estimate can be derived for the components of the partition of 



s'/., so that there holds 



e" 

E4"|Af -Afc|=0(l).|AQ(ze)|. (5.10) 



Therefore, relying on the inductive ssumption, from (|5.9p - (l5.10p we recover the 
desired estimate p.l5p . which completes the proof of the proposition. □ 



6 Conclusion 

Here we briefly describe how to get the proof of Theorem 11.21 following the ideas 
contained in |10j and relying on the results established in the previous section. 
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step 1. 

We use the partition of waves of an approximate solution u' into 

primary waves { , A } , 

secondary waves {Vk^^k] ■> 



provided by Proposition 15 . II to construct a piecewise constant approximation w = 
w(t,x) of u'^{t,x) in a time interval [me,ne] that enjoys the following properties 
(see [TOl Section 4]). 

1. The wave fronts in w are of two kinds, primary and secondary. 

2. Each primary front originates at t = me and ends at t = ne; 

3. There is a one-to-one correspondence between primary fronts and primary 
waves {Vk}- In particular, the primary front corresponding to ykin^jj) has 
constant size 's'j^{m,j) and, in view of Proposition 15. 11 joins with a segment 
the points {me,je) and {ne^l(^n,j,k,h)£) of the {t,x) plane. 

4. The left and right states of the primary front corresponding to y^{m,j), say 
u'l'^{m,j), u'^'^'{m, j), are always related by 

Uk^'imJ) = Tk[u'l'^ (m, j)] (s^(m,j)) . 
Moreover, there holds 

w{me) = u'^(me) . 

5. Let u^{t), and u^{t) be the left and right state of a secondary front X/3(t) 
of w at time t e [m£,ne]. Then, letting denote the set of all pairs of 
crossing primary waves in it^ (i.e. all pair of waves y^{m,j), y^'/(?Ti,j') for 
which j < j\k > k' and t{n,j,kji) > ^{n,j' ,k' ,h')), there holds 



J2\uf{t)-uf,{t)\^Oil) 




],k,h CW 



= 0(1) • jAT"'"! , 

where the summand on the left hand side runs over all secondary fronts in 
w(t), while the second summand on the right hand side runs over all pairs 
of crossing primary waves in m^. 

6. All secondary fronts travel with speed 2, strictly larger than all characteristic 
speeds. 
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step 2. 



Using the same arguments of [TUl Section 5], relying on p.l2p . p.l7p . p.l4p . 
p.lSp . one can prove that 



\S(n-,n}ewime) - w{ne)\\^^^ 



= 0(1)- lAT^ 



1 + log(n — to) 



n — m 



{n-m)e, (6.1) 



\u^ne) ~ w{ne)\^^ = 0(1) • |AT"'"| • [n - m)e 



where S'(„_m)eit'(w£) is the semigroup trajectory related to (|l.ip with initial datum 
wlrne) = u^{nne) evaluated at time t = (n — m)e. 

Step 3. 

Now, as in |TD1, Section 6], let T = me + e', for some to S N, < e' < e, and fix 
a positive constant p > 2e. Them, we define inductively integers = toq < TOi < 
. . . < TO^ = TO in this way. Assuming given, then 

1. if T(TOi£) — T((TOi + 1)£) < p, let TOi+i be the largest integer less or equal 
to m such that (mi+i — TOi)£ < p and T(mie) — T(TOi+i£) < p; 

2. if T(mi£) — T((TOi + 1)£) > p, set m.i+i = rui + 1. 

On every interval [to^e, m^+iE] where 1. holds, we construct a piecewise constant 
approximation of u'^ according to Step 1, and using (|6.ip we derive 



|M^(m,+i£) - S(^„ii+i-mi)eU%m^e)\\^, 

-m^m. + l] , 1 + l0g(w,+ i - TO^) 



0(1) 



AT''' 



( TOi+i — mi)e . 



(6.2) 



On the other hand, on each interval [TOi£,TOi+i£] where 2. is verified, by the 
Lipschitz continuity of and applying (|1.17p we find 



[^^(to^+ie) - S'(,„^_^^_,„^)5U^(to,£)||j^i = 0(1) • £ . 



(6.3) 



Hence, observing that the cardinality of both classes of intervals 1.-2. is bounded 
by 0(1) • from ^^-^^ we finally deduce 



\\u'{T)-Stu\\^,^0{1) 
which yields (|1.13p choosing p = ^/e ■ log | log£|. 



p+-log-+£ l + - 
p £ \ p 



□ 
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Remark 6.1 By the same observations of Remark \3. 61 one deduces that the con- 
clusion of Theorem l 1 . S\ remains valid if we assume that the flux function F is C^'^ 
and that, for each k-th NGNL characteristic family, the linearly degenerate man- 
ifold Mk in (|1.8p is the union of a finite number of connected manifolds Mk,h, 
that are either N—\- dimensional as in (H), or N -dimensional and in this case the 
following conditions hold. The flux function F is onQ,\Mk, the vector field r^ 
is transversal to the boundary of Aik,h, and letting d^A4k,h, d~A4k,h denote the 
connected components of the boundary of M.k,h where rk points towards fl\Mk.h 
and Mk^h7 respectively, the one-sided second derivatives of Xk (see Remark \3.6\) 
satisfy 

V+(VAfc •rfc)(w) TfcCu) < yued+Mk, 

V-(VAfe-rfe)(u)-rfc(M) <0 ^ued-Mu, 

or the opposite inequalities. Indeed, if we again add the term 2 ^ p,g p | in 

the interaction potential Qq defined by (|4.7p - (|4.8p . for every wave Sa containing 
several contact discontinuities j^, . . . , the estimates stated in Lemma \4-6\ 
continue to hold, provided that 

inf{s>0 : Rk[u]{s) ed+Mk,h, ued-MkM} >0, (6.5) 

which is certainly true up to a possible slight restriction of the domain Q. Relying 
on Lemma \4-6\ one then establishes Proposition \5.1\ with the same arguments of 
Section\^ and thus conclude as above. 
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